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Abstract. The DAMA/Libra experiment has recently confirmed the annual modulation signal 
obtained in the earlier DAMA/Nal experiment, providing strong evidence that they have actually 
detected dark matter. We examine the implications of this experiment for the mirror dark matter 
candidate. We show that mirror dark matter successfully explains the latest DAMA/Libra data, 
whilst also being consistent with the null results of other direct detection experiments. 
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I'm interested in talking about some recent developments in efforts aimed at the direct 
detection of dark matter. There are a number of on-going experiments, but there is 
one experiment which is particularly interesting because it seems that they've actually 
detected dark matter. This is the DAMA/Libra experimentlll] and its predecessor the 
DAMA/Nal experiment [1]. 

These DAMA experiments eliminate the background by using the annual modulation 
signature. The idea[3] is very simple. The interaction rate must vary periodically since it 
depends on the Earth's velocity, v#, which modulates due to the Earth's motion around 
the Sun. That is, 

R(v E ) =R(v Q ) + (j^J Av e cos co(t -to) (1) 

where Ave — 15 km/s, co = 2n/T (T = 1 year) and to = 152.5 days (from astronomical 
data). The phase and period are both predicted! This gives a strong systematic check 
on their results. Such an annual modulation was found[l] at the 8.2c Confidence level, 
with r, to measured to be: 

T = 0.998 ±0.003 year 

t = 144 ±8 day (2) 

Clearly, both the period and phase are consistent with the theoretical expectations of halo 
dark matter. These are strong reasons to believe that the DAMA people have detected 
dark matter. The data, together with the cosine prediction is given in figure 1. 

What type of theory could explain these results? As a first try, one could take a 
simple minded approach to the question of dark matter. One could imagine that the 
laws governing dark matter could be identical to the laws governing ordinary matter. By 
this I mean that the dark matter could belong to a sector which is an exact duplicate of 




the ordinary matter sector, so that the Lagrangian is: 

J? = J?SM(e,u,d,y 1 ...)+J?sM(e',u' 1 d\y',...) (3) 

Such a theory can also be motivated from a symmetry reason if left and right handed 
chiral fields are interchanged in the extra sector. In this way space-time parity symmetry 
can be realized as a symmetry of nature, and for this reason we call the particles in the 
extra sector mirror particles. The standard model extended with a mirror sector was first 
studied in ref.[4] and shown to be a consistent extension of the standard model (for a 
review and more complete list of references see ref.[5|]). 

Observe that P' ,e' ,He' etc would have masses identical to their ordinary matter 
counterparts. They would also be stable and thus potentially dark matter candidates. In 
order to make this type of theory consistent with the successful big bang nucleosynthesis 
(BBN) measurements and successful large scale structure (LSS), we need to assume that 
in the early Universe, i.e. when the temperature was around 1 MeV and less, that T' < T. 

In fact, some studies suggest that we require T'/T ~ 0.3 Jg]. 

The idea as it stands doesn't explain any dark matter direct detection experiment 
since the ordinary and mirror particles interact with each other only via gravity. A 
relevant question is: can ordinary and mirror particles interact with each other non- 
gravitationally? That is, can we add any interaction terms consistent with renormaliz- 
ability and the symmetries of the theory? The answer is YES - but only two terms are 
possible[|4|]: 
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FIGURE 2. Rutherford elastic scattering of mirror nuclei off ordinary nuclei induced by photon - mirror 
photon kinetic mixing of strength e. 

where F^ v (i^ v ) is the ordinary (mirror) U(l) gauge boson field strength tensor and 
(0') is the electroweak Higgs (mirror Higgs) field. Both of these terms can lead to 
interactions between the ordinary and mirror particles. It turns out that the coupling of 
the Higgs to protons and neutrons is too weak for the Higgs mirror Higgs quartic term to 
be important in dark matter experiments. This leaves us with the photon mirror photon 
kinetic mixing term, which it turns out, can do the job fzLS]. In fact this term will lead 
to Rutherford-type elastic scattering (figure 2) of ordinary nuclei off mirror nuclei (for 
the relevant cross-section, form factors, and other technical details, see ref.[7|,|8D). 

To make contact with the direct detection experiments we need to know something 
about the chemical composition and distribution of particles in the galactic halo. 

Chemical composition: 

H' \He' are expected to be produced in the early Universe, and in fact we expect a larger 
He' mass fraction, Y^ e i > Yfj e , ifT'<T. Elements heavier than He', such as 0\ Ne' .... 
are expected to be produced in mirror stars. 

Distribution: 

To explain the rotation curves in spiral galaxies, we know that the dark matter needs to be 
roughly spherically distributed in spiral galaxies. Given the upper limit on compact star 

sized objects in the halo from MACHO searchers (roughly f macho ~ 0.2 — 0.3 depending 
on the assumptions), we then expect the mirror particles to be distributed predominately 
as a gaseous spherical halo surrounding the collapsed disk of ordinary matter[9(]. A 
dissipative dark matter candidate like mirror matter can only survive in an extended 
spherical distribution without collapsing if there is a substantial heating mechanism to 
replace the energy lost due to radiative cooling. It turns out, that there is such a heating 
mechanism, which I will comment on, at the end of this talk. 

Anyway, the gas of mirror particles will have a Maxwellian velocity distribution, 



= e 2 




(5) 



where the index i labels the particle type [i = e' ,P',He',0' , ...]. In the non-interacting 



one species WIMP case, T = jinv^ ot , (where v rot ~ 250 km/s is the rotational velocity in 
the Milky Way). In the mirror matter case where we have several particle species with 
significant self interactions, we must examine the condition for hydrostatic equilibrium, 
which balances the central gravitational attraction with the pressure gradient: 

£ = -p« (6) 



dr 



where 



P = J>/r, p = £m,n ; , g =°J r pdV = V -^. (7) 

Here, n,,p, T are the number density of species i (i = e',P',He', ...), mass density and 
Temperature, which is assumed to be isothermal for simplicity. It is straightforward to 
solve this equation, leading tod]: 

T = \mv 2 rot (8) 

where m = Y* m i n i/Y* n i ls the average mass of the particles in the halo. Thus we see that 
the velocity dispersion of the particles in the dark matter halo depends on the particular 
particle species and satisfies: 



Vo['1=vL- (9) 
nii 

Note that if m\ in, then Vq[i\ <C vjl ot . That is, heavier mirror particles will have a 
very narrow velocity distribution, which turns out to be a key feature in explaining why 
DAMA sees a signal and the other experiments do not (within this framework). In fact, 
given that the kinetic energy of halo mirror particles in the Earth's reference frame is 
E = \mv 1 rot ~ keV for m,- = m# e , it follows that with DAMA's 2 keV recoil energy 
threshold that they are indeed only sensitive to elements heavier than He, that is the sub 
dominant component, naively expected to be mirror oxygen O' , and from the above, will 
have narrow velocity distribution Vq[0'] <C v^ ot . 

Fixing the normalization of the annual modulation gives a measurement of e, which 
turns out to be e ~ 10 9 [7]. A value for e of around this magnitude is consistent 
with experimental and astrophysical constraints (for a review, see. ref. niQTI ). and also 
consistent with early Universe cosmology bounds (successful BBN and LSS)lll|]. 

Importantly, the DAMA/Libra experiment provides new information because they 
have enough data to begin to examine the dependence of their signal on the recoil 
energy of the target nuclei, which is something that they measure. We define the annual 
modulation amplitude S m by: 

dR(v E ) dR(v & ) fddR/dE R \ 

= + S m cos 0)(t- 1 ) . (10) 
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FIGURE 3. Energy dependence of the cosine modulation amplitude, S m , for four illustrative cases: 
A' = Si',v rot = 170 km/s (solid line) A' = Mg' ,v m = 195 km/s (long-dashed line), A' = Ne',v rot = 230 
km/s (short-dashed line), A' = O' . v wt = 280 km/s (dotted line). In each case e is fixed so that the mean 
amplitude is 0.0129 cpd/kg/keVee. Also shown are the DAMA/Nal & DAM A/LIBRA combined data 
from figure 9 of ref . OJJ] . This figure assumes that the mass of the halo is dominated by He'. 

How do we expect S m to vary with Er! For large Er 3> \mv 2 rot , only particles in the 
tail of the Maxwellian distribution could produce such a large Er scattering event in 
the Earth based detector. On June 2nd, when ve, is a maximum, we expect more large 
Er scattering events, since there will be a lot more mirror particles which are energetic 
enough (in Earth's reference frame) to produce a given large Er scattering. Thus, in the 
large Er region, we expect S m to be positive and tending to zero as Er — » °o. At low 
Er <C jinVy 0t , all the particles in the halo are energetic enough to produce a scattering 
event, but since do IcIEr 1/v 2 , we expect less events when ve is a maximum, so that 
S m < at low Er. To summarize, we expect S m to have an Er dependence which goes 
to as Er — > °°, rises to a peak as Er moves from infinity towards the body of the 
distribution, and changes sign at low Er. 

Note that the position of the peak, which is essentially given by kinematics, will give a 
measure of the mass of the dark matter particle. In ref. [7] I have done the analysis, taking 
into account all the required things, such as energy resolution of the detector, quenching 
factors etc. Assuming a halo dominated in mass by H' /He' with an A 1 subcomponent 
(where A' ~ O' is the expectation), I find a predicted shape (see figure 3) for S m which 
agrees nicely with the shape measured by DAMA/Libra. This is a significant test of the 
theory. 

Allowing for a range of possible v rot (most recent astronomical data suggest lfl2ll v rot = 
254 km/s), leaves a band of allowed parameter space, shown in figure 4. Assuming v rot = 
254 ± 30 km/s, a galactic halo consisting of a dominant H'.He' with A' subcomponent 
can fit the data for A' with mass 18 ±4 GeV, which is compatible with the naive 
expectation of O' (given that oxygen is the most abundant ordinary element after H/He). 
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FIGURE 4. DAMA/Libra allowed region together with the 95% C.L. exclusion limits from the 
XENON10[13], CDMS/Ge and CDMS/Si I0] experiments (the regions above the exclusion contours 
are the disfavoured region). 



It turns out that ALL of this parameter space is consistent with the null results of the 
other experiments. The 95% C.L. limits of the most sensitive of these null experiments, 
are also indicated. 

These results are in sharp contrast to popular models involving weakly interacting 
massive particles (WIMPs). Elastic scattering of standard WIMPs give a very poor 
fit to the data ni5h . The basic reasons as to why elastic scattering with mirror dark 
matter works, and elastic scattering of standard WIMPs doesn't has to do with their 
basic differences: a) the velocity distribution of O' dark matter is much narrower than 
for WIMPs, Vq[0'] < v 2 rot , while for standard WIMPs Vq = v 2 rot , and b) Rutherford 
scattering has a cross section do /cIEr 1 jE\, while for standard WIMPs do /<1Er is 
Er independent (excepting the Er dependence on the form factors). I would also argue 
that mirror dark matter is simpler and more elegant than models with standard WIMPs, 
so in a sense its not surprising that experiments have come out in favour of mirror dark 
matter, but of course I might be biased! 

This mirror matter interpretation of the dark matter detection experiments will be 
tested further by DAMA - as they collect more data their statistical error will reduce. 
Furthermore the DAMA people plan to upgrade their experiment replacing their PMTs 
with the aim of lowering their Er threshold. This will be particularly useful, as they 
should see the change in sign of their modulation S m , predicted in figure 3. Experiments 
like XENONlotfH] and CDMSlH] have E R threshold too high to see the same mirror 
dark matter component which DAMA is detecting. They could still find a positive signal 
if there is a heavier component, such as an Fe' component. Such a component should be 
there at some level. 

The inferred value of £ ~ 10~ 9 seems quite interesting for supernova physics. In 



particular for such values of £ about half of the total energy emitted in Supernova 
explosions will be in the form of light mirror particles (v' e „ T , e f± , /) Hill . This implies 
a heating of the halo (principally due to the e f± component), of around: 

LhZat-in ~ \ x 3 x 10 53 oyJ^ ~ 10 44 erg/s for Milky Way (11) 

It turns out that this matches (to within uncertainties) the energy lost from the halo due 
to radiative cooling |@]: 

Lh ener g y-out = A J r n 2 e AKr 2 dr ~ 10 44 erg/s for Milky Way. (12) 

In other words, a gaseous mirror particle halo can survive without collapsing because 
the energy lost due to dissipative interactions is replaced by the energy from ordinary 
supernova explosions. Presumably there is some detailed dynamical reasons maintaining 
this balance, which of course, may be difficult to elucidate due to the complexity of that 
particular problem. 
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